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Amendments to the Claims: 

This listing of claims will replace all prior versions, and listings, of claims in the application: 
Listing of Claims: 



that binds to a protein when coordinate data of the protein of single chain or plural chains is 
given, the apparatus comprising: 



structural change in consideration of dynamic behavior using an induced-fit parameter 
reflecting induced fit on the coordinate data of the protein and selects post-structural-change 
protein coordinate data; 

a spatial point designating unit that designates a spatial point at which 
superposition with the ligand is to be conducted[[,]] from the post-structural-change protein 
coordinate data selected by the post-structural-change protein coordinate data selecting unit; 

an interaction function calculating unit that calculates an interaction function 
when the protein and the ligand bind to each other using the spatial point designated by the 
spatial point designating unit and a ligand coordinate data of the ligand; and 

a ligand evaluating unit that evaluates the ligand that binds to the protein based 
on the interaction function calculated by the interaction function calculating unit. 

2. (Currently Amended) The ligand screening apparatus according to claim 1, 
wherein the interaction function calculating unit calculates the interaction function using 
Scor e (ij) Sscore(U) shown in Formula 1 [[.]]: 



1 . (Currently Amended) 



A ligand screening apparatus which screens for a ligand 



a^post-structural-change protein coordinate data selecting unit that effects 



when i is not equal to j 




v 



[Formula 1] 



when i is equal to j 



ax(l-/?) 
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[[(]]wherein dy S is a distance between i-th spatial point and j-th spatial point in the 
target protein[[.]] i dij C is an interatomic distance between i-th atom and j-th atom in [[the]] a 
compound[[.]] a a is a coefficient for making Sscore(i j) the maximum value when the group 
of spatial points in the target protein and the compound completely overlap with each 
other[[.]] a p is a coefficient for giving a threshold value by which it can be defined as 
"overlapping"[[)]]. 

3. (Currently Amended) The ligand screening apparatus according to claim 1, 
wherein the interaction function calculating unit further comprises an interaction function 
optimizing unit that carries out optimization so as to make the score of interaction function 
maximum. 

4. (Currently Amended) The ligand screening apparatus according to claim 3, 
wherein the interaction function calculating unit further comprises: 

an interaction energy optimizing unit that calculates interaction energy with 
the protein for the superposed ligand after optimization of the interaction function by the 
interaction function optimizing unit, and optimizes the interaction energy while finely 
adjusting conformation of ligand 3D structure data. 

5. (Currently Amended) The ligand screening apparatus according to claim 4, 
wherein the ligand evaluating unit further comprises: 

a reevaluating unit that execut e s causes execution of the interaction function 
calculating unit after largely changing conformation of ligand 3D structure data following 
optimization by the interaction energy optimizing unit, and reevaluates the ligand that binds 
to the protein based on the interaction function calculated by the interaction function 
calculating unit. 

6. (Currently Amended) The ligand screening apparatus according to claim 1, 
wherein in calculation of any one of the induced-fit parameter and the post-structural-change 
protein coordinate data or both, the post-structural-change protein coordinate data selecting 
unit calculates a normal mode for the protein coordinate data, determines intensity of 
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fluctuation of each amino acid, and conduct conducts a molecular dynamic calculation using 
the intensity of fluctuation as a constraint condition. 

7. (Currently Amended) The ligand screening apparatus according to claim 6, 
wherein the post-structural-change protein coordinate data selecting unit calculates a 
fluctuation value of a_dihedral angle of a main chain according to normal mode calculation, 
and conducts a molecular dynamic calculation while setting the fluctuation value as a 
coefficient of force K in the molecular dynamic calculation shown by Formula 2 or Formula 
3[[.]]1 

Erot = Krot (<J) - <j)0) 2 [Formula 2] 

[[(]]wherein Erot represents energy of a dihedral angle of a main chain atom in a 3D 
structure of a protein^.]],, § represents a dihedral angle of the main chain atom[[.]] a $0 
represents astandard value of the dihedral angle of the main chain atom[[.]] i [[Here,]] when a 
value of Krot (a coefficient of force ) is large, <j> is constrained by <|>0[[.)]]. 

Epos = Kpos (r - r0) 2 [Formula 3] 

[[(]]wherein Epos represents position energy of the main chain atom in a 3D structure 
of a protein[[.]] A r represents a_coordinate of the main chain atom[[.]] a rO represents a standard 
value of the coordinate of the main chain atomff.]^ [[Here,]] when a value of Kpos £a 
coefficient of force) is large, r is constrained by r0[[.)]]. 

8. (Withdrawn) The ligand screening apparatus according to claim 1, wherein 
the interaction function calculating unit uses the interaction function to which a dynamic 
property function representing dynamic property of protein is added as "elastic energy". 

9. (Currently Amended) The ligand screening apparatus according to claim 8, 
wherein the interaction function calculating unit adapts "U collision" as "elastic energy" 
which is a function shown by Formula 4 in consideration of local flexibility of the protein[[.]] 
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M N 

u collision = Z2M'>>) 

<p (/, y ) = Kcollision * (Rcollision (i 9 j)-R) 2 

[Formula 4] 

[[(]]wherein M is anumber of atoms in an active site that prohibit collision, N is a 
number of atoms of the ligandr|".l], [[When]] interatomic distance R between an i-th atom of a 
main chain or a side chain with a little dynamic behavior in an active site, and a j-th atom in 
the ligand is not more than collision distance "Rcollision(ij)", (|)(i,j) is calculated[[)]]. 

10. (Currently Amended) The ligand screening apparatus according to claim 1, 
wherein the interaction function calculating unit uses the interaction function to which a 
normal mode analysis result or secondary structure determination result of the protein is 
added as a dynamic property function that represents a dynamic property of the protein. 

1 1 . (Currently Amended) A ligand screening method which screens for a ligand 
that binds to a protein when coordinate data of the protein of single chain or plural chains is 
given, the method comprising: 

apost-structural-change protein coordinate data selecting step that effects,, 
with a computer, structural change in consideration of dynamic behavior using an induced- fit 
parameter reflecting induced fit on the coordinate data of the protein and selects post- 
structural-change protein coordinate data; 

a spatial point designating step that designates a spatial point at which 
superposition with the ligand is to be conducted[[,]] from the post-structural-change protein 
coordinate data selected by the post-structural-change protein coordinate data selecting step; 

an interaction function calculating step that calculates an interaction function 
when the protein and the ligand bind to each other using the spatial point designated by the 
spatial point designating step and a ligand coordinate data of the ligand; and 
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ajigand evaluating step that evaluates the ligand that binds to the protein 
based on the interaction function calculated by the interaction function calculating step; 

wherein the interaction function calculating step calculates the interaction function 
using Sscore(ij) shown in Formula 1: 



Sscore («,./) = £ 

u 



wherein d \ f is a distance between i-th spatial point and i-th spatial point in the target 
protein, djj C is an interatomic distance between i-th atom and i-th atom in a compound, a is a 
coefficient for making Sscore(i r i) the maximum value when the group of spatial points in the 
target protein and the compound completely overlap with each other, ft is a coefficient for 
giving a threshold value by which it can be defined as "overlapping" . 

12. (Cancelled) 

13. (Currently Amended) The ligand screening method according to claim 1 1, 
wherein the interaction function calculating step further comprises an interaction function 
optimizing step that carries out optimization so as to make the score of interaction function 
maximum. 

14. (Currently Amended) The ligand screening method according to claim 13, 
wherein the interaction function calculating step further comprises: 

an interaction energy optimizing step that calculates interaction energy with 
the protein for the superposed ligand after optimization of the interaction function by the 
interaction function optimizing step, and optimizes the interaction energy while finely 
adjusting conformation of ligand 3D structure data. 
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ax[exp{-(^-^) 2 }-/? 



[Formula 1] 



when i is equal to j 
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15. (Currently Amended) The ligand screening method according to claim 14, 
wherein the ligand evaluating step further comprises: 



calculating step after largely changing conformation of ligand 3D structure data following 
optimization by the interaction energy optimizing step, and reevaluates the ligand that binds 
to the protein based on the interaction function calculated by the interaction function 
calculating step. 

16. (Currently Amended) The ligand screening method according to claim 11, 
wherein in calculation of any one of the induced-fit parameter and the post-structural-change 
protein coordinate data or both, the post-structural-change protein coordinate data selecting 
step calculates a normal mode for the protein coordinate data, determines intensity of 
fluctuation of each amino acid, and conduct conducts a molecular dynamic calculation using 
the intensity of fluctuation as a constraint condition. 

17. (Currently Amended) The ligand screening method according to claim 16, 
wherein the post-structural-change protein coordinate data selecting step calculates a 
fluctuation value of a^dihedral angle of amain chain according to normal mode calculation, 
and conducts a molecular dynamic calculation while setting the fluctuation value as a 
coefficient of force K in the molecular dynamic calculation shown by Formula 2 or Formula 



[[(]]wherein Erot represents energy of a dihedral angle of a main chain atom in a 3D 
structure of a protein[[.]] a (() represents a dihedral angle of the main chain atom[[.]] A <j)0 
represents a_standard value of the dihedral angle of the main chain atom[[.]] a [[Here,]] when a 
value of Krot (a coefficient of force) is large, § is constrained by <j>0[[.)]L 



^reevaluating step that executes causes execution of the interaction function 



3[[.]]: 




[Formula 2] 




[Formula 3] 
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[[(]]wherein Epos represents position energy of the main chain atom in a 3D structure 
of a protein[[.]] r represents ^coordinate of the main chain atom[[.]] a rO represents a standard 
value of the coordinate of the main chain atom[[.]] i [[Here,]] when a value of Kpos £a 
coefficient of force) is large, r is constrained by rO[[.)]L 

1 8. (Withdrawn) The ligand screening method according to claim 1 1 , wherein 
the interaction function calculating step uses the interaction function to which a dynamic 
property function representing dynamic property of protein is added as "elastic energy". 

19. (Currently Amended) The ligand screening method according to claim 18, 
wherein the interaction function calculating step adapts "U collision" as "elastic energy" 
which is a function shown by Formula 4 in consideration of local flexibility of protein[[.]]; 

M N 

U co |ii s i on = ^j^ J ( P 0 J ) rT? t A i 

m j~\ [Formula 4] 

<p (i 9 j) = Kcollision * (Rcollision (i, j) - R) 2 

[[(]]wherein M is a number of atoms in an active site that prohibit collision, N is a 
number of atoms of ligand[[.]L [[When]] interatomic distance R between an i-th atom of a 
main chain or a side chain with a little dynamic behavior in an active site, and a j-th atom in 
the ligand is not more than collision distance "Rcollision(ij)", c()(i,j) is calculated^)]],. 

20. (Currently Amended) The ligand screening method according to claim 1 1 , 
wherein the interaction function calculating step uses the interaction function to which a 
normal mode analysis result or secondary structure determination result of the protein is 
added as a dynamic property function that represents a dynamic property of the protein. 

2 1 . (Currently Amended) A non-transitory computer-readable medium having a 
program which makes a computer execute a ligand screening method which screens for a 
ligand that binds to a protein when coordinate data of the protein of single chain or plural 
chains is given, the method comprising: 

a post-structural-change protein coordinate data selecting step that effects 
structural change in consideration of dynamic behavior using an induced-fit parameter 
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reflecting induced fit on the coordinate data of the protein and selects post-structural-change 
protein coordinate data; 

a spatial point designating step that designates a spatial point at which 
superposition with the ligand is to be conducted, from the post-structural-change protein 
coordinate data selected by the post-structural-change protein coordinate data selecting step; 

an interaction function calculating step that calculates an interaction function 
when the protein and the ligand bind to each other using the spatial point designated by the 
spatial point designating step and a ligand coordinate data of the ligand; and 

ajigand evaluating step that evaluates the ligand that binds to the protein 
based on the interaction function calculated by the interaction function calculating step; 

wherein the interaction function calculating step calculates the interaction function 
using Sscore(i r i) shown in Formula 1 

Sscore(ij) = ^ 

if 



wherein d \ f is a distance between i-th spatial point and i-th spatial point in the target 
protein, djj C is an interatomic distance between i-th atom and i-th atom in a compound, a is a 
coefficient for making SscoreH j) the maximum value when the group of spatial points in the 
target protein and the compound completely overlap with each other, ft is a coefficient for 
giving a threshold value by which it can be defined as "overlapping" . 

22. (Cancelled) 

23. (Currently Amended) The non-transitory computer-readable medium having a 
program according to claim 21, wherein the interaction function calculating step further 
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when i is equal to j 
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comprises an interaction function optimizing step that carries out optimization so as to make 
the score of interaction function maximum. 

24. (Currently Amended) The non-transitory computer-readable medium having a 
program according to claim 23, wherein the interaction function calculating step further 
comprises: 

an interaction energy optimizing step that calculates interaction energy with 
the protein for the superposed ligand after optimization of the interaction function by the 
interaction function optimizing step, and optimizes the interaction energy while finely 
adjusting conformation of ligand 3D structure data. 

25. (Currently Amended) The non-transitory computer-readable medium having a 
program according to claim 24, wherein the ligand evaluating step further comprises: 

a reevaluating step that e x e cutes causes execution of the interaction function 
calculating step after largely changing conformation of ligand 3D structure data following 
optimization by the interaction energy optimizing step, and reevaluates the ligand that binds 
to the protein based on the interaction function calculated by the interaction function 
calculating step. 

26. (Currently Amended) The non-transitory computer-readable medium having a 
program according to claim 21, wherein in calculation of any one of the induced-fit parameter 
and the post-structural-change protein coordinate data or both, the post-structural-change 
protein coordinate data selecting step calculates a normal mode for the protein coordinate 
data, determines intensity of fluctuation of each amino acid, and conduct conducts a 
molecular dynamic calculation using the intensity of fluctuation as a constraint condition. 

27. (Currently Amended) The non-transitory computer-readable medium having a 
program according to claim 26, wherein the post-structural-change protein coordinate data 
selecting step calculates a fluctuation value of a dihedral angle of a main chain according to 
the normal mode calculation, and conducts the molecular dynamic calculation while setting 
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the fluctuation value as a coefficient of force K in the molecular dynamic calculation shown 
by Formula 2 or Formula 3 [[.]]: 

Erot = Krot (<|> - <j>0) 2 [Formula 2] 

[[(]]wherein Erot represents energy of a dihedral angle of a main chain atom in a 3D 
structure of a protein[[.]] i <() represents the dihedral angle of the main chain atom[[.]] a $0 
represents a_standard value of the dihedral angle of the main chain atom[[.]]* [[Here,]] when 
a value of Krot (a coefficient of force) is large, is constrained by <t>0[[.)]L. 

Epos = Kpos (r - rO) 2 [Formula 3] 

[[(]] wherein Epos represents a position energy of the main chain atom in 3D structure 
of a protein[[.]] a r represents a coordinate of the main chain atom[[.]] a r0 represents a standard 
value of the coordinate of the main chain atom[[.]] a [[Here,]] when a value of Kpos £a 
coefficient of force) is large, r is constrained by rO[[.)]L 

28. (Withdrawn) The program according to claim 21, wherein the interaction 
function calculating step uses the interaction function to which a dynamic property function 
representing dynamic property of protein is added as "elastic energy". 

29. (Currently Amended) The non-transitory computer-readable medium having a 
program according to claim 28, wherein the interaction function calculating step adapts "U 
collision" as "elastic energy" which is a function shown by Formula 4 in consideration of 
local flexibility of protein[[.]]; 

M N 

u collision = Z !>('>./) 

1=1 7=1 

<p (i, j) = Kcollision * (Rcollision (i, j)-R) 2 

[Formula 4] 

[[(]]wherein M is ajiumber of atoms in an active site that prohibit collision, N is a 
number of atoms of ligand[[.]] a [[When]] interatomic distance R between an i-th atom of a 
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main chain or a side chain with a little dynamic behavior in active site, and a j-th atom in the 
ligand is not more than collision distance "Rcollision(ij)", (|)(i J) is calculated[[)]]. 

30. (Currently Amended) The non-transitory computer-readable medium having a 
program according to claim 21, wherein the interaction function calculating step uses the 
interaction function to which a normal mode analysis result or secondary structure 
determination result of the protein is added as a dynamic property function that represents a 
dynamic property of the protein. 

31. (Cancelled) 
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